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ABSTRACT 


This  report  summarizes  one  year  of  theoretical  and 
applied  research  on  propagation  of  seismic  waves  and  tech¬ 
niques  for  analyzing  data.  The  main  objectives  were  to 
determine  che  frequency  and  energy  of  seismic  signatures, 
and  investigate  attenuation,  patterns  of  azimuthal  radia¬ 
tion  from  source  regions,  and  methods  of  determining  the 
type  of  motion  at  the  source.  Natural  and  artificial 
sources  were  studied  to  develop  diagnostic  aids  t‘or  dis¬ 
tinguishing  between  earthquake*  and  underground  nuclear 
detonations.  Equipment  for  selection,  reformatting,  and 
digital-to-analog  conversion  for  digitally  recorded  LASA 
data  was  constructed  and  is  being  checked  out.  Several 
approaches  for  using  the  parallel  computational  capabili¬ 
ties  of  optics  for  LASA  data  were  developed.  A  study  of 
background  noise  and  reciprocity  for  teleseismic  events 
as  recorded  on  the  bottom  of  a  large  fresh  water  lake  has 
commenced  with  the  emplacement  of  three-component  seismo¬ 
meters  in  Lake  Superior.  Array  data  have  been  useo  for 
crustal  studies  on  the  Eastern  United  States.  Digital 
mode  filtering  was  investigated.  A  perturbation  theory 
for  seismic  sources  was  developed. 
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COLLECTION  AND  ANALYSIS  OF  SEISMIC  WAVE 
PROPAGATION  DATA 


1 

INTRODUCTION 

A  program  of  research  into  the  collection  and  analyses 
of  data  concerning  the  propagation  of  seismic  waves  was  con¬ 
ducted  at  the  Geophysics  Laboratory  of  Willow  Run  Laboratories 
a  unit  of  the  University  of  Michigan's  Institute  of  Science 
and  Technology,  during  the  one-year  period  from  1  June  1966 
through  31  May  1967.  This  research  was  supported  by  the 
Advanced  Research  Projects  Agency  and  was  monitored  by  the 
Air  Force  Office  of  Scientific  Research  under  Contract  AF49 
(638)-17S9.  It  was  a  continuation  of  studies  begun  under 
sponsorship  of  Contracts  AF49(638)-91  1  ,  AF49(638)-1078,  and 
AF49(638)-1170, 

The  main  objectives  of  the  program  were  to  determine  the 
frequency  and  energy  of  seismic  signatures,  study  attenuation, 
investigate  the  patterns  of  azimuthal  radiation  from  source  re¬ 
gions,  and  study  methods  of  determining  the  type  of  motion  at 
the  source.  Natural  and  ai'^ificial  sources  were  studied  in  or¬ 
der  to  develop  diagnostic  aids  to  distinguish  between  earth¬ 
quakes  and  underground  nuclear  detonations.  Equipment  for  se¬ 
lection,  reformatting,  and  digi tal - to-analog  conversion  for  the 
large-aperture  seismic  array  (LASA)  data  has  been  constructed 
and  is  being  checked  out.  Several  approaches  for  using  the 
parallel  computational  capabilities  cf  optics  were  studied  for 
their  potential  contributions  to  immediate  identification  for 


digital  processing  of  seismic  disturbances.  An  investigation 
of  seismic  background  noise  and  the  reciprocity  of  signal  am¬ 
plitudes  for  teleseismic  events  as  recorded  on  the  bottom  of 
a  large  fresh  water  lake  has  commenced  with  the  emplacement 
of  ocean  bottom  three-component  seismometers  in  Lake  Superior. 
Array  data  from  the  Michigan  Basin  through  North  Carolina 
have  been  used  for  crustal  studies.  Jigital  computation  of 
mode  filtering  to  separate  recordings  of  body  waves  from  Ray¬ 
leigh  waves  has  been  investigated.  A  perturbation  theory  for 
seismic  sources  has  been  developed. 
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RESULTS  AND  DISCUSSIONS 

2.1.  DESCRIPTION  OF  CONVERTOR 

The  large  aperture  seismic  array  (LASA)  data  convertor, 
designed  to  select  and  reformat  recordings  from  the  LASA  digi¬ 
tal  magnetic  tapes  is  completely  constructed.  However,  the 
convertor  is  not  yet  in  use  because  final  optimization  of  tape 
readout  remains  to  be  accomplished.  The  logic  circuits  have 
been  checked  and  debugged.  Originally  designed  to  select  LASA 
data  for  the  inv'stigation  of  optical  processing,  the  convertor, 
within  reasonable  limits,  is  adaptable  for  input  to  any  type 
of  digital  or  analog  data  handling  equipment.  The  output  of 
the  convertor  can  be  recorded  on  digital  or  analog  magnetic  tape, 
or  punched  on  paper  tape,  or  presented  on  a  CRT.  During  the 
recording  process  continuous  monitoring  can  be  performed.  While 
channels  or  subarray  sums  are  selectively  recorded  and  changed 
at  will,  presentations  of  the  output,  the  record  number,  parity 
checks,  etc.,  are  given  visually,  and  the  input  tape  can  be  ad¬ 
vanced  or  reversed  to  select  desired  i^ecorded  intervals. 

Both  the  projected  optical  investigations  and  the  logical 
design  of  the  LASA  data  convertor  a  e  described  in  the  final 
report  (5178-64-F)  on  the  previous  contract.  Two  new  possible 
approaches  to  the  optical  investigations  have  been  developed 
since.  The  first,  convolution-function  spatial  filtering,  is 
described  in  a  publication  in  Applied  Optics  (see  appendix  I); 
the  second,  developed  primarily  under  another  contract,  is  a 
simple  optical  method  of  obtaining  sectional  autocorrel ograms 
and  retrocorrelograms ,  which  recently  have  been  shown  to  be 
interpretively  beneficial  in  oil  exploration  seismology  [1]. 
Further  work  during  this  period  has  also  shown  that  lensless 
correlography  can  be  suitably  performed  with  the  input  as  a 
CRT  pattern. 
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Since  the  above-mentioned  report  described  the  logical 
scheme  of  the  LASA  data  convertor,  the  operational  character¬ 
istics,  described  by  the  designer,  Rowland  H.  McLaughlin, 
fol low. 

2.1.1.  LASA  DATA  CONVERTOR 

LASA  is  a  long-range  reismic  detection  system,  consisting 
basically  of  a  large  number  of  seismic  detectors  arratiged  in 
geometric  patterns.  By  proper  correlation  of  the  data  * rom 
the  individual  detectors,  the  array  may  be  so  steered  as  to 
have  maximal  sensitivity  in  any  direction.  Proper  correlation 
of  the  data  will  also  achieve  large  improvements  in  signal-to- 
noise  ratio  over  that  available  from  individual  detectors. 

One  of  the  outputs  available  from  LASA  will  consist  of  magnetic 
tape-recordings  of  the  outputs  from  the  individual  detectors, 
in  digital  form.  The  recording  format  provides  that  a  sample 
of  the  data  from  each  detector  will  be  taken  every  50  ms  and 
the  digitized  output  from  the  detectors  will  be  multiplexed 
together  with  status  information  and  date-time  code'  into  a 
form  compatible  with  IBM  digital  magnetic  tape-recordings.  The 
main  facets  of  the  format  pertinent  to  the  following  discussion 
are ; 

(1)  One  sample  of  data  from  an  individual  seismometer 
is  called  a  data  word*  and  consists  of  18  bits**. 

The  data  word  for  a  specific  seismometer  occupies 
a  unique  position  in  each  frame. 

(2)  A  frame  is  one  sample  of  data  from  each  seismome¬ 
ter  plus  weather  and  telemetry  information.  The 
telemetry  orovides  information  for  calibration 
and  maintenance. 


*A  word  is  an  ordered  set  of  characters  that  has  at  least 
one  meaning,  and  is  stored  or  transferred  by  the  processor  as 
a  unit. 

**A  bit  is  the  smallest  unit  of  information  recognized  by 
the  processor. 
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(3)  A  record  will  include  two  frames  plus  one  header, 
and  occasionally  one  trailer. 

(4)  A  header  containing  auxiliary  data,  including 
time,  date,  record  identification,  and  "trailer 
present"  indicator  initiates  each  record. 

(5)  The  trailer  consists  of  auxiliary  data  to  be  re¬ 
corded  with  an  occasional  record,  and  includes  the 
status  of  each  detector,  the  status  of  each  sub¬ 
array,  and  similar  information. 

The  LASA  data  convertor  (figs.  1  and  d)  has  been  constructed 
for  the  specific  purpose  of  reformatting  LASA  magnetic  d;,ta 
tapes.  It  is  simpler  and  faster  than  a  general  p.-pose  computer 
primarily  because  only  thj  desired  data  is  read  from  the  tape. 

The  unit  contain--  a  digita  i-to-analog  convertor  as  an  integral 
component,  and  hence  the  data  extracted  from  the  master  tape 
may  be  provided  to  the  user  in  either  digital  or  analog  form.  A 
block  diagram  of  the  system  for  extracting  the  appropriate  data 
from  the  LASA  recordings  is  shown  in  figure  3. 

The  unit  shown  in  figure  3  performs  three  basic  functions. 
First,  appropriate  data  on  the  LASA  magnetic  tape  is  recognized 
by  its  position  in  the  tape-recording  format  and  transferred 
to  a  core  memory  unit;  second,  a  signal  is  generated  by  the 
daca-transfer  unit  which  may  be  used  for  auxiliary  functions 
later,  such  as  placing  the  LASA  data  in  the  proper  position  on 
the  multiple-channel  film  (this  infor.mation  is  stored  in  the 
core  memory  along  with  the  LASA  data);  and  third,  the  LASA  data 
and  its  auxiliary  information  are  read  from  the  core  memory  on 
command.  An  internal  command  is  generated  at  50  ms  intervals 
(to  restore  the  real-time  correlation  of  the  data  samples) 
whenever  stored  data  are  converted  to  cnalog  signals. 
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FIGURE  1.  I  ASA  Da td-Convertor . 

power  supply  rack. 
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FIGURE  2.  LASA  Da td-Ccnver tor . 


Data  processing  rack. 
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Auxiliary  function  of  the  data-transfer  unit  are: 

(1)  Recognizing,  counting,  and  displaying  the  pro¬ 
gression  of  records  during  the  reproduction  of 
a  magnetic  tape 

(2)  Controlling  the  magnetic  tape-transport  motion 
either  automatically  by  use  of  decoded  signals 
from  the  tape  cr  by  manual  control 

(3)  Displaying  selected  data  contained  in  header, 
trailer,  telemeter  codes,  or  data  words 

(4)  Decoding  and  converting  the  time  code  to  a  form 
suitable  for  recording  on  film 

Figure  4  shows  the  controls  used  to  select  data  which  is 
to  be  read  from  the  master  tape  into  the  core  memory  of  the 
coiivertor.  As  mentioned  previously  data  samples  from  each 
seismometer  occupy  the  same  unique  position  in  each  frame. 

The  selection  controls  are  thumbwheel  cwitches  arranged  in 
groups  of  three.  Each  thumbwheel  is  individually  adjustable  so 
that  each  selector  group  may  be  rapidly  set  for  any  number  from 
0' 0  to  999.  Fifty  selector  groups  are  provided  and  are  numbered 
sequentially  and,  hence,  fifty  data  samples  may  be  selected  by 
dialing  the  selectors  to  correspond  with  the  seismometer  number 
whose  output  is  desired.  Using  a  particular  selector  group, 
say  number  five,  also  generates  a  code  unique  to  the  selector 
group  chosen,  and  this  code  is  stored  in  the  magnetic  core 
memory  along  with  the  seismometer  sample.  This  auxiliary  fea¬ 
ture  is  useful  in  certain  analog  applications  such  as  position¬ 
ing  the  converted  data  sample  on  a  film  recording  or  CRT  display. 
Two  additional  switches  are  associated  with  each  of  the  fifty 
data  selector  groups.  The  data-bias  switch  permits  a  bias  level 
to  be  generated  for  drawing  spacing  lines  on  a  film  recording. 
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FIGURE  4.  Infornation  sel ecti on . panel 


The  frame  switch  al.ows  the  selection  of  data  from  frame  1, 
frame  ?.  ^  or  both  depending  on  the  setting  chosen. 

Consider,  for  example,  that  it  is  desired  to  select  the 
delta  in  both  frames  generated  by  the  center  seismometer  in 
subarray  D4,  Consulting  the  LASA  magnetic  tape  format,  D4  is 
the  thirteenth  subarray  in  the  sampling  sequence,  and  the  cen¬ 
ter  seismometer  in  each  subarray  is  sampled  first.  This  means 
that  an  appropriate  selector  group,  say  number  10,  is  set  to 
013.  The  data-bias  switch  is  set  to  "data"  and  the  frame 
switch  is  set  to  "both".  This  setting  's  shown  in  figure  4. 

A  similar  procedure  is  followed  for  each  data  sample  desired. 
The  result  of  such  a  setting  is  that  eacn  time  the  m.=»ster  tpae 
is  incremented  one  record,  two  samples  of  data  from  the  chosen 
seismometer  are  stored  in  a  known  location  in  the  core  memory, 
and  are  available  for  readout  into  any  desired  device. 

The  present  design  of  the  convertor  reads  the  data  out  of 
the  core  at  50  ms  ’ntervals  and,  by  means  of  the  integral  digi¬ 
tal  -  to-anal  og  convertor,  converts  the  data  to  analog  form. 

Other  types  of  readout  ir  analog  form,  for  example,  into  an 
analog  emputer  or  into  ac  analog  magnetic  recorder,  can  also 
be  made.  Present  facilities  only  perm  *:  time  division  multi¬ 
ple)'  readout  of  multiple  samples,  however.  Of  course,  a  single 
sample  may  be  selected  and  a  continuous  readout  obtained.  The 
digi tal-to-analog  converter  used  has  provision  for  demultiplex¬ 
ing  six  channels  of  data  into  six  separate  outputs,  and  if  this 
feature  is  desired  the  demultiplexing  hardware  coulo  be  pur¬ 
chase  and  installed. 

If  digital  output  is  desired,  such  as  punched  paper  tape 
or  digital  magnetic  tape,  additional  equipment  could  be  added 
to  provide  this  type  of  readout. 


In  addition  to  the  fifty  data-word  selectors,  four  addi- 
tion.n  selectors  are  connected  t.o  storage  register-display 
combinations.  These  four  units  shown  in  figure  5  each  contain 
18  bits,  and  each  bit  controls  a  lamp.  These  registers  are  pri¬ 
marily  provided  for  reading  out  the  header  information.  How¬ 
ever,  since  the  selectors  associated  with  these  registers  are 
also  adjustable,  they  may  be  set  so  as  to  read  out  and  display 
additional  data,  such  as  the  output  from  any  specific  seismo¬ 
meter,  or  the  contents  of  various  portions  of  the  trailer. 

The  number-four  display  register  is  provided  with  one 
additional  option.  A  permanently  wired  weather  decoder  may  be 
connected  to  this  display.  This  provision  is  not  merely  a  con¬ 
venience  but  a  requirement,  since  the  weather  is  sub-multi¬ 
plexed  across  four  records.  The  details  of  the  multiplexing 
are  contained  in  the  LASA  Magnetic  Tape  Format. 

The  number-one,  number-two,  and  number*  three  display  re¬ 
gisters  are  also  equipped  with  adjustable  digital  comparators 
that  permit  the  operator  to  search  for  and  control  additiTnal 
functions  from  particular  information. 

A(i  auxiliary  feature  of  the  LASA  data  convertor  that  is 
provided  for  operator  convenience  is  the  two-way  BCD-binary- 
number  translator  (fig.  6).  This  device  allows  the  operator 
to  insert  a  BCD  number,  for  example,  and  read  the  binary  equi¬ 
valent  from  the  display.  Thus,  a  means  of  rapid  conversion 
from  one  number  system  to  the  other  is  provided.  This  unit 
may  also  be  used  to  read  out  and  display,  for  operator  interpre¬ 
tation,  the  information  btored  in  the  core  memory,, 

A  record  counter  may  be  seen  in  figure  5.  This  counter 
counts  the  number  of  records  processed  since  the  beginning  of 
the  tape.  This  unit  also  contains  provisions  for  presetting 
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FIGU^^E  6.  Number  translator  panel 


the  counter  to  any  arbitrary  number.  Associated  with  the 
record  counter  are  two  selectors  that  can  provide  control  of 
the  tape's  motion.  The  selectors  may  be  set  to  any  number 
contained  in  the  record  counter.  Two  such  selectors  a'^e  pro¬ 
vided,  which  permits  the  unit  to  be  set  so  as  to  shuttle  be¬ 
tween  two  preset  numbers.  In  addition,  manual  starting, 
manual  forward-backward,  and  manual  rewinding  controls  are 
provided  (fig.  7). 

2.2.  COHERENT  OPTICAL  PROCESSING 

The  investigation  of  coherent  optical  processing  per  se 
had  been  concluded  in  the  previous  contract.  However,  to  ob¬ 
tain  visual  output  of  the  fine  frequency  structure,  operational 
use  of  the  processor  was  employed  on  recordings  from  three 
events:  LONGSHOT,  EARLY  RISE,  and  GREELEY. 

As  described  in  the  section  on  earthquake  data,  the  pro¬ 
blem  on  LONGSHOT  was  to  determine  the  nature  of  a  "pre-PcP" 
arrival,  that  is,  whether  the  arrival  was  a  differently  re¬ 
fracted  P,  or  was  a  deep  reflection  similar  to  the  PcP.  To 
make  this  interpretation  Fourier  transforms  of  six  second  in¬ 
tervals  of  the  "pre-PcP"  arrival,  the  P  and  the  PcP  were  made. 

The  frequency  distribution  of  the  "pre-PcP"  was  compared  to 
those  of  the  P  and  PcP.  The  fine  frequency  structure  available 
from  the  optical  Fourier  transform  was  instrumental  in  indi¬ 
cating  that  the  pre-PcP  arrival  was  most  propably  a  deep  reflec¬ 
tion. 

Portions  of  the  EARLY  RISE  data  were  processed  to  aid  in 
comparison  of  the  consistency  of  the  fine  structure  from  sta- 
tion-to-station  across  the  Michigan  Basin,  and,  along  with  the 
GREELEY  data,  for  a  comparison  with  the  results  of  mode  filtering 
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FIGURE  7.  Tape  transport  control  panel. 


and  to  determine  the  effects  of  crustal  variation  on  the  fine 
frequency  structure. 

A  new  method  of  optical  spatial  filtering — correlation 
function  filtering — was  developed.  In  this  technique  the  two- 
dimensional  cross-  or  auto-correlation  of  two  functions  is 
operated  upon  by  ‘‘elective  attenuation.  One  of  the  functions 
is  simultaneously  imaged,  resulting  in  a  "time-domain"  equi¬ 
valent  of  complex  frequency  filtering.  It  is  planned  to  inves¬ 
tigate  the  utility  of  this  technique  for  processing  LASA  data. 
The  technique  is  described  in  an  article  in  Applied  Optics  (see 
appendi x  I ) . 

2.3.  MODE  FILTERING 

Work  has  been  done  on  the  technique  of  distinguishing  P 
from  Rayl ei gh-type  wave  motion  using  filtering  techniques  per¬ 
taining  to  the  four  quadrants  of  the  Fouri er-transf orm  plane. 
Originally  done  optically  on  synthetic  and  empirical  data,  the 
mode-filtering  process  has  been  adapted  to  digital  computation. 
The  digital  approach  allows  g»'eater  flexibility  in  the  inter¬ 
pretation  and  graphic  presentation  of  the  filtered  output. 

Using  a  program  which  generates  synthetic  seismic  input,  fil¬ 
ters  it,  and  reconstitutes  the  filtered  output,  contour  maps 
of  signal-to-noise  ratios  for  various  values  of  the  geometric 
and  signal  parameters  were  produced. 

Mathematical  analysis  has  produced  analytic  expressions  for 
M(a,t)  and  N(a,t)  (c.f.,  Tifteenth  Quarterly  Progress  Report  on 
contract  AF49(638)-1 1 70) .  The  expressions  are  sufficiently  com¬ 
plex  to  preclude  the  feasibility  of  numeric  -.rialysis  by  hand. 

A  computer  program  was  written  to  calculate  vai‘vies  of  M(a,t)  and 
N(a,t)  so  that  they  can  be  presented  in  tabular  or  a  variety  of 
graphical  forms. 
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Hand  computation  is  successful  in  one  case,  however, 
namely  that  in  which  signal  and  noise  are  of  the  same  fre¬ 
quency  and  cc-d i recti  or  a  1 ,  and  in  which  the  carrier  is  rotated 
so  ti.at  Sp  =  Sp  =  Tr/2.  Let  the  input  vertical  signal-to-noise 
ratio  be  c"/c,  the  iriput  radial  signal-to-noise  ratio  be  d'/d; 
tlien  the  output  vertical  signal-to-noise  ratio  of  the  signal 
filtered  with  0(6, oj)  (c.f.,  fifteenth  Quarterly  Progress  Report) 
i  s 


fs 
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Z  out 


\JY.9b 

— Lji _  J 

2.95  , 


When  c  =  c"  =  3/  2  and  d-'d'  =  l. 


which  is  in  close  agreement  with  the  experimentally  obtained 
value  of  about  3  (c.f.,,  Se  ,th  Semiannual  Technical  Summary 
Report  on  contract  AF49(638)-1  1  70) . 

2.4.  THEORETICAL  STUDIES  IN  ELASTIC  WAVE  SCATTERING  (I.  Mclvor) 

A  perturbation  method  ha'  oeen  developed  for  obtaining  the 
scattered  solution  due  to  imperfections  on  a  plane  boundary. 

It  has  applications  to  a  number  of  wave  propagation  problems  of 
seismic  interest.  The  basic  geometry  is  shown  in  figure  8.  The 
plane  y=0  is  the  nominally  plane  boundary  between  two  elastic  re¬ 
gions.  The  imperfections  in  the  boundary  are  defined  by  g=g(x), 
whicli  represents  the  deviation  of  the  actual  boundary  from  the 
pi ane  y=0. 
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Figure  8 


The  problem  presently  being  investigated  is  the  reflection 
of  a  plane  compress! onal  wave  by  a  free  surface  (i.e.,  I'egion  2 
has  zero  rigidity).  A  formal  solution  has  been  obtained  for  the 
t.ee  surface  motion  for  small  arbitrary  deviations  g.  When 
e.  ’uated  in  specific  cases,  the  results  may  be  of  interest  in 
several  applications;  e.g.,  the  effect  of  imperfections  on  the 
apparent  angle  of  eme^’gence  can  be  determined. 


The  same  approach  will  give  results  for  a  number  of  simi¬ 
lar  problems.  A  solution  can  be  obtained  for  the  effect  of  im¬ 
perfections  on  surface  waves.  This  is  of  interest  in  noise 
studies  and  possibly  woulo  have  significance  in  the  attenuation 
of  surface  waves.  Results  obtained  by  the  method  for  the  cwo- 
layered  region  shown  in  the  figure  could  be  applied  to  reflec¬ 
tions  from  the  Mohorovicic  discontinuity.  The  effect  of  imper¬ 
fections  in  the  "Moho"  on  such  reflections  may  contribute  to 
the  attenuation  of  reflect-^d  wave  phases.  It  is  anticipated 
that  these  extensions  of  the  problem  presently  being  considered 
will  be  investigated  in  the  near  future. 

To  illustrate  the  basic  ideas  of  the  pe'"turbation  method, 
its  formulation  in  the  context  of  the  free  surface  reflection 
problem  is  outlined  here.  Except  for  the  statement  of  the 
boundary  conditions,  the  formulation  is  identical  for  the  other 
problems  referred  to  in  the  previous  pa-agraph. 

The  displacement  components  associated  with  steady-state 
motion  in  the  x-y  plans  can  be  expressed  in  terms  of  two  poten¬ 
tial  functions.  This  representation  is  with  the  time  factor 
e’“^  suppressed 


3x  3y 

3y  3x 

The  dilatational  and  shear  potentials  satisfy 

V'4)*  +  =  0 

=  0 


(1) 


(2; 
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wi  th 


ki"  =  oiVci^.  ka"  =  ujVc2^  (3) 

where  Ci  and  C2  are  respectively  the  dilational  and  shear  wave 
speeds  given  by 

Ci^  =  (A+2y)/p,  C2’  =  u/p  (4) 

in  which  X  and  y  are  Lame's  constants  and  p  is  the  mass  density. 

The  solution  of  (2)  representing  the  reflection  of  a  plane 
compres s i ona 1  wave  by  a  free  plane  boundary  is  given  in  any 
standard  text  on  elastic  wave  propagation.  A  quantity  associated 
v/ith  this  solution  is  denoted  here  with  a  subscript  zero,  e.g., 
♦q.  For  the  present  problem  we  introduce  the  scattered  poten¬ 
tials  as  the  difference  between  the  actual  potentials  and 
Thu? 

a  4,  +  4,^,  »  ij)  +  Vp  (5) 


where  (|)  and  4)  are  the  scattered  potentials.  The  present  problem 
is  then  defined  by  the  boundary  value  problem 


7^4,  +  k;H  =  0 
7^4^  +  k2^4^  =  9 


(6) 


where  the  solutions  must  satisfy  the  radiation  condition  and  on 
the  boundary  y=g{x)  the  stress  condi^-ions 


on+a  n=-(a  n+a.  n) 
XX  xy  y  '  xo  X  xyo  y' 


0  n+an=-(o  n+o  n) 
xy  X  y  y  '  xyo  x  yo  y' 


(7) 
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in  which  n  ,  n  are  the  dire.cion  cosines  of  the  outer  normal 
X  y 

to  the  boundary.  We  let  g(x)  =  af(x)  where  a  denotes  a  char¬ 
acteristic  amplitude  of  the  i moerf ecti on .  Thus  f(x)  may  be 
taken  as  0(1). 

The  major  difficulty  of  the  problem  is  that  the  boundary 
is  not  a  coordinate  curve.  This  can  be  avoided  by  an  appro¬ 
priate  transformation  at  the  expense  of  complicating  the 
governing  differential  equations.  If  the  imperfections  are 
small,  however,  the  transformed  equations  will  involve  a  small 
parameter.  Thus,-  they  are  amenable  to  perturbation  techniques. 

The  wavelength  of  the  incident  wave  is 

K  =  ki/271  (8) 

We  restrict  ourselves  to  small  imperfections,  i.e.,  a/K  <<1,  and 
introduce  the  small  parameter 

e  =  2ira/K  (9) 

With  this  a  convenient  coordinate  transformation  is 

5  =  2Trx,/K,  n  =  2iTy/K  -  cc(c)  (10) 

where  c(5)  ■  f(K5/2Tr).  In  the  coordinates  the  free  boundary 
is  given  by  q-O.  Introducing  the  transformation  (10)  into  equa¬ 
tions  (6)  yields 

4>"  +  (l+c^^'M*!)"  -  c(2^'<p''  +  ?"<f)  +  4>  =  0 

(11) 

<p"  +  (1+c2;'2)4;"  -  c(2r-i>'’  +  +  a^4>  =  0 

where  the  prime  represents  d/d^  and  the  dot  8/8q,  and 

=  Cl 

The  boundary  conditions  are  still  give.)  by  (7),  but  are  new  to  be 
evaluated  at  r,=0. 
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We  next  introduce  the  expansions 

(})  =  C({)i  +  e^!})2  +  ... 

ip  =  eip  I  +  + 


eo  +  e ‘ o  +  .  .  . 

Xi  X2 


yi  ■  y2 

eo  +  c^a  + .  ,  . 
xyi  xy2 

u/ ( 2~/ K)  =  GUi  +  e^U2+... 
v/(2Tr/K)  =  eVi  +  e^V2  +  ... 


(12) 


The  direction  cosines  are  functions  of  e.  Their  power  series 
expansions  are 

n^  =  e;'(l  -  1/2(ec')'+...)  (13) 

ny  =  -1  +  l/2(er)2  +  ... 

The  stresses  o^  etc.,  associated  with  the  perfectly  smooth 
boundary  will  a?so  be  functions  of  e  when  evaluated  at  ri*0. 

They  cat.  be  expanded  in  power  series  to  yield 


-  0 


-a  +  eo*  +e^o„^  +. 

X  x„ 

0  0  0 

=  eo*  +  e^a  ^  , 

^0  ^0 


xy, 


=  eo*  +e^o 

xyo  xy^+, 


(14) 


where,  of  course,  the  coefficients  on  the  right  hand  side  are 
known  functions. 
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When  equations  (12),  (13),  and  (14)  are  introduced  into 
(11)  and  (7),  we  obtain  a  series  of  boundary  value  problems 
giving  successively  higher  order  approximations.  The  various 
order  displacement  and  stress  quantities  can  be  expressed  in 
terms  of  the  associated  potential  functions  by  introducing 
equations  (12)  into  the  displacement  representation  and  the 
stress-strain  relations.  We  record  here  only  the  first  order 
approximation.  It  is 


<i>i"  +  4>i  +  4>i  =  0 
'Pi"  +  Pi  +  =  0 


(15) 


with  boundary  conditions  at  ri=0 


where 


and 

with 


xyi 


yi 


a 


^yo 


=  pki  (-Y<Jii  +  2<pi"  +  2p:''} 
X  1 

Oy^  =  ukiM-Y4>i  +  24)V  -  2iJ;r*) 

=  yki  +  'j^V  -  'j^T') 

xy  1 

Ui  =  I  '  +  i|;i,  V 1  =  1  ■  'Pi 


Y  =  A/U  • 


(16) 


(17) 


(18) 


A  solution  to  this  first  order  approximation  has  been  ob¬ 
tained  for  arbitrary  functions  c;.  It  has  the  form  of  convolution 
integrals  of  the  imperfection  t;  with  furrctions  of  ^  and  n.  These 


24 


functions  are  expressed  as  complex  integrals.  Although  alge¬ 
braically  involved,  they  can  be  explicitly  integrated  by  means 
of  contour  integration  and  asymptotic  expansions.  They  are 
presently  being  evaluated.  When  this  is  carried  out,  results 
for  specific  problems  can  be  obtained.  These  will  be  given 
in  future  reports. 

2.5.  FIELD  MEASUREMENTS 
2.5.1.  PROJECT  EARLY  RISE 

Eight  University  of  Michigan  field  teams  participated  in 
the  recording  of  the  Lake  Superior  “Early  Rise"  shots  during 
July  1966.  One  hundred  and  twenty-one  sites  were  located  along 
a  southeast  line  from  the  Keweenaw  Peninsula,  Upper  Michigan  to 
the  North  Carolina  Coastal  Plains.  These  stations  partially  re- 
ve»'sed  a  seismic  refraction  profile  obtained  from  the  ECOOE  and 
CHASE  shots.  The  station  spacing  from  a  <  800  km  was  approxi¬ 
mately  10-15  km.  At  greater  distances  a  station  spacing  of 
approximately  25  km  was  utilized.  At  the  larger  distances  two 
to  eight  shots  or  more  were  recorded  at  the  same  site  in  an 
effort  to  obtain  reliable  data.  This  was  necessitated  because 
of  high  seismic  background  noise  with  resulting  poor  signal-to- 
noise  ratios.  Matched  three-component  short-period  instruments 
with  FM  magnetic  tape  recorders  were  used  at  each  site. 

To  date,  of  the  308  recordings  obtained,  296  have  been 
processed  and  detailed  frequency  analyses  made  of  146  recordings. 
A  preliminary  report  of  the  EARLY  RISE  data  containing  travel 
times,  distances,  amplitudes,  and  predominant  frequencies  was 
prepared  and  submitted  to  the  sponsor  previously.  Iso-particle 
velocity  maps  have  been  prepared  for  the  vertical  and  longi¬ 
tudinal  components  of  first  P,  maximum  P  and  maximum  shear- 
surface  wave  arrivals  at  each  site.  Detailed  attenuation  studies 
are  currently  being  completed  for  these  data.  Array  processing 
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techniques  were  employed  at  five  sites  using  the  multiple  shots 
recorded  at  the  same  site.  Impressive  signal  - to-noise  improve¬ 
ments  (a  factor  of  four)  were  obtained  at  the  most  cistant  site 
(1680  km)  using  this  technique.  However,  uncertainties  in  the 
preliminary  shot  point  locat'^ons  degraded  this  technique  so 
work  was  suspended  pending  receipt  of  corrected  shot  locations. 
Now  that  final  shot  point  locations  are  available  this  investi¬ 
gation  v.'ill  be  continued. 

Travel  time  studies  of  the  EARLY  RISE  'ecordings  disclosed 
an  intermediate  crustal  layer  with  a  velocity  of  7.7  km/sec.  A 
mantle  velocity  of  8.3  km/sec  and  an  upper  crustal  velocity  of 
6.2  km/sec  were  also  indicated.  Scattered  travel  time  resi¬ 
duals  across  the  Appalachian  Mountains  make  the  determination 
of  the  crustal  thickness  uncertain.  Lower  travel  time  resi¬ 
duals  on  the  western  flank  of  this  range  may  indicate  a  change 
in  n?;  le  velocity  (8.5  km/sec)  or  may  indicate  a  thinning  of 
the  crustal  layer.  The  larger  residuals  observed  across  the 
crest  of  the  range  may  be  due  to  a  mountain  root  system  such  as 
that  observed  for  the  Rocky  Mountains. 

2.5.2.  ECOOE  AND  CHASE  SERIES 

I  so- par t i cl e- vel oci ty  maps  have  been  prepared  for  the  ECOOE 
shots  recorded  during  1965.  Six  stations  recorded  the  CHASE  VII 
event.  Detailed  frequency  analyses  have  been  completed  for 
these  recordings.  A  composite  report  is  under  preparation  which 
will  present  the  results  of  the  attenuation  and  propagation 
studies  using  data  of  the  EARLY  RISE,  ECOOE,  CHASE,  1962  North 
Carolina,  and  1953-1964  Lake  Superior  shots.  A  presentation  of 
this  report  is  scheduled  for  the  lUGG  meetings  at  Zurich 
(September-October  1967). 

2.5.3.  NTS  EVENTS 

Fifteen  NTS  events  were  re'orded  in  Michigan  during  the 
past  year.  Six  of  these  events  were  recorded  at  2  to  5  stations 
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while  the  remainder  were  recorded  at  only  the  Botanica’’  Gardens 
well  station.  Frequency  analyses  were  performed  for  9  NTS 
events  . 

2.5.4.  EARTHQUAKE  DATA 

Numerous  earthquakes  have  been  recorded  on  the  intermediate 
depth  well  system  (Botanical  Gardens)  during  the  past  year  using 
the  anticipator  signal  threshold  triggering  device  that  has  been 
discussed  in  previous  reports.  The  ESSA  Preliminary  Determina¬ 
tion  of  Epicenter  (PDE)  cards  received  during  this  period  have 
been  scanned  and  a  catalogue  of  the  listed  events  recorded  on 
the  well  system  was  prepared.  Unfortunately  time  and  funds  have 
not  permitted  analyzing  much  of  this  data.  One  particular  task 
which  is  under  study  is  the  determination  of  the  spectral  con¬ 
tent  of  compressional  waves  generated  by  deep  focus  earthquakes. 
It  is  anticipated  that  tnis  study  will  aid  in  the  interpreta¬ 
tion  of  the  crustal  data  (effect  of  local  crustal  structure  on 
spectral  content  of  P)  obtained  from  the  underwater  HE  shots 
detonated  during  the  last  fi\e  years. 

A  pt"eliminary  examination  of  passbandi  from  the  well  sys¬ 
tem  recording  of  a  deep-focus  earthquake  (2/15/67,  Brazil-Peru 
border,  Mag.  6.2,  depth  600  km)  revealed  significant  high- 
frequency  (above  6  cps)  energy  In  the  P,  PeP  and  SeP  phases. 
Epicentral  distance  was  53°. 

The  spectral  content  of  local  and  regional  earthquakes 
is  also  under  study.  Records  of  the  Columbus,  Ohio  earttiquake 
of  April  8,  1967  and  the  Buffalo,  New  York  quake  of  June  13, 

1967  are  being  studied.  The  spectra  of  these  earthquakes  will 
be  compared  with  that  of  an  event  recorded  prior  to  the  CHASE 
VII  shot  in  1966.  This  event  was  well  recorded  at  several  of 
the  mobile  stations  deployed  for  the  CHASE  shot.  Additional 
records  have  been  obtained  from  other  stations  in  the  region, 
the  epicenter  has  been  located  along  the  Ohio-Indiana  border 
at  2  distance  of  about  200  km  from  Ann  Arbor.  The  magnitude 


of  this  event  is  about  1.0  on  the  Richter  scale  and  the  large 
amplitude  of  S  indicates  that  it  may  be  an  earthquake.  If  so, 
it  is  possible  that  many  of  the  unexplained  and  often  unnoticed 
small  high-frequency  arrivals  found  on  short-period  records 
from  this  region  may  be  mi  croearthqua!,es .  The  high  background 
levels  make  higher  station  density  impc.ative  if  these  events 
are  to  be  adequately  identified. 

An  effort  to  compare  the  P  and  PcP  spectra  of  the  LONGSHOT 
event  as  recorded  at  stations  in  North  Michigan  led  to  the  dis¬ 
covery  of  an  arrival  about  15  seconds  ahead  of  PcP  and  travel¬ 
ing  with  the  same  apparent  velocity.  A  search  is  being  con¬ 
ducted  to  identify  this  phase  on  other  records  and  to  determine 
whether  it  represents  a  deep  refracted  path  or  a  reflection 
from  just  above  the  core-mantle  boundary.  Records  of  the  LONG- 
SHOT  ev^nt  and  of  intermediate  and  deep  earthquakes  in  the  dis¬ 
tance  range  55®  to  65®  have  been  obtained  from  the  VELA  Seismo- 
logical  Center  for  this  purpose. 

2.5.5.  LAKE  SUPERIOR  LAKE  BOTTOM  SEISMOMETER  STUDY 

In  April  1967  two  ocean  bottom  seismometer  packages  were 
obtained  on  loan  from  the  University  of  California  for  the  pur¬ 
pose  of  determining  the  seismic  background  noise  on  the  bottom 
of  a  large  fresh  water  lake  and  to  determine  the  reciprocity  of 
signal  amplitudes  for  teleseismic  events.  These  units  had  pre¬ 
viously  been  used  in  the  Pacific  Ocean  for  gathering  seismic 
data  .'or  short  periods  of  time. 

Earlier  in  the  year  personnel  of  the  Geophysics  Laboratory 
made  a  trip  to  the  Institute  of  Geophysics  and  Planetary  Physics, 
University  of  California  to  coordinate  the  details  of  the  in¬ 
strumentation  and  underwater  techniques.  In  March  1967,  Mr.  John 
B'aumler  spent  two  weeks  at  LaJolla  for  the  purpose  of  becoming 
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intimately  acquainted  with  these  units  as  they  had  been  used 
in  the  ocean  and  the  methods  of  dispensing  and  recovering 
them  at  sea. 

As  received,  these  units  were  apparently  ready  for  use 
under  the  same  conditions  which  the  University  of  California 
utilized  them  in  the  ocean.  For  the  University  of  Michigan's 
pregram  in  Lake  Superior  it  was  necessary  to  make  some  c^'anges 
to  encompass  the  frequency  passband  desired  and  the  length  of 
time  they  must  operate  ''nderwater  unattended.  In  addition  to 
the  above,  additional  facilities  were  installed  to  simplify 
data  reduction. 

With  the  first  of  several  drops  in  Lake  Superior  using 
these  packages  due  in  the  month  of  June  1967  it  appeared  at 
the  time  that  there  would  be  no  problems  in  implementing  this 
schedule.  However,  this  did  not  p  ove  to  be  the  case. 

The  fi^st  modification  made  was  in  the  internal  clock  un^t 
(SEIKO)  used  to  put  time  markers  on  the  magnetic  recording  tape. 
It  was  discovered  that  random  clock  pulses  were  being  generated 
and  recorded  on  tape.  In  order  to  eliminate  this  problem  it 
was  necessary  to  incorporate  a  filter  circuit  in  the  clock  and 
also  in  the  tape  recorder  capstan  motor  in  order  to  prevent 
the  commutation  and  speed  regulating  circuit  from  interfering 
with  the  pulse:  produced  by  the  clock  circuit.  The  clocks  were 
also  modi  fifed  to  provide  both  minute  and  hour  markers  in  addi¬ 
tion  to  the  original  10  cps  and  1  cps  pulses.  The  original 
10  cps  and  1  cps  signals  were  rectangular  pulses  approximately 
20  milliseconds  wide.  The  "relay  box"  circuitry  of  the  clock 
was  modified  such  that  the  10  cps  signal  was  converted  into  a 
square  wave.  In  addition,  the  minute  and  hour  pulses  were 
brought  out  in  an  "or"  configuration  instead  of  a  "coincidci.ee" 
configuration.  Both  the  minute  and  hour  pulses  were  used  to 
turn  "on"  a  transistor  "gate"  circuit  into  which  was  also  fed 
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the  10  cps  signal.  When  put  on  one  channel  of  the  tane  re¬ 
corder,  there  then  appeared  a  one-minute  time  marker,  one- 
second  in  width  filled  with  the  10  cps  square  wave.  When  the 
hour  marker  appeared,  its  duration  was  one-minute  filled  with 
the  10  cps  square  wave.  This  was  done  to  give  a  more  accurate 
determination  of  time.  With  the  above  additions,  search  and 
’oqging  time  have  been  greatly  reduced  in  data  re_overy. 

The  underwater  seismograph  units  were  equipped  with  Texas 
Instrument's  multichannel  type  I  parametric  seismic  amplifiers. 
Input  attenuators  for  these  amplifiers  were  built  and  installed 
with  the  consent  of  University  of  California  personnel.  The 
parametric  amollTier  has  a  basic  gain  of  60  db  followed  by  a 
high-cut  filter.  Two  additional  fixed  20  db  stages  of  ampli- 
fi  ion  are  also  included.  This  facilitates  a  triple  gain 
feature  for  each  seismometer  of  60,  80  and  100  db.  These  three 
outputs  arc  recorded  on  individual  channels  on  the  tape  recorder. 

As  used  by  the  University  of  California,  the  high-cut  fil¬ 
ters  were  set  for  a  10-second  period  with  a  24  db/octave  cutoff. 

A  much  higher  cutoff  frequency  was  necessary  for  the  Lake 
Superior  program.  Hence,  values  recommended  ^'or  a  10  cps  filter 
were  installed.  The  besv  that  could  be  done  with  these  values 
was  a  slope  of  about  16  db/octave.  After  considerable  experi¬ 
mentation  it  was  found  that  the  radial  change  in  resistance 
value  of  the  RC  time  constants  in  the  filter  circuit  affected 
the  'ransistor  biasing  and  feedback  to  sucn  a  degi^ee  that  the 
circuit  did  not  operate  to  within  the  stated  specifications. 

By  modifying  the  circuit  values  in  the  high-cut  filter,  tne  de¬ 
sired  24  db/octave  cutoff  at  17  cps  was  obtained  with  a  phase 
shift  between  channels  of  no  more  than  5  degrees  overall.  On 
the  apolication  of  a  signa’  to  tne  input  of  the  amplifier  it 
was  found  that  as  che  amplitude  of  the  signal  was  increased  and 
the  tLree  outputs  of  the  amplifier  were  monitored  (60,  80,  and 
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100  db),  the  first  output  to  overload  and  distort  was  the  100 
db  position.  As  the  signal  input  war;  increased  by  about  6  db 
the  60  db  outpu;.  became  distorted  and  consequently  this  in 
turn  immediately  distorted  the  80  Jb  output.  It  was  found  that 
the  output  of  the  high-cut  filter  was  being  loaded  by  the  10,000 
ohm  input  impedance  of  the  tape  recorder  bias  oscillator. 

A  transistor  emitter  follower  circuit  was  added  between  Ihe 
high-cut  filter  output  and  the  inout  of  th  "low"  channel  on 
the  tape  recorder.  This  not  only  prevented  loading  of  the  high- 
cut  filter  but  also  provided  greater  dynamic  range  to  the  ampl'- 
f  i  er . 

The  three  Teledyne  Earth  Sciences  Division  Model  SD213 
"Ranger"  type  seismometers  contained  in  each  of  the  ocean  bot¬ 
tom  units  were  apparently  no  longer  "tuned"  to  the  two-second 
period  as  s’own  in  the  original  specifications  sheets.  These 
seismometers  are  "detuned"  by  the  proximity  of  magnetic  material 
and/or  their  relation  to  each  other  when  placed  close  together. 
It  was  necessary  to  "charge"  each  of  thf  seismometers  to  a 
longer  period  than  desired,  then  "discharge"  them  individually 
until  they  were  at  the  correct  period  so  that  they  could  be 
phase  matched  to  within  a  few  degrees.  The  orthogonal  set  of 
seismometers  in  the  first  unit  was  tuned  with  a  natural  period 
of  1.1  seconds.  In  addition  to  the  tuning  of  the  seismometers 
by  "charging"  them  to  frequency,  a  6  yFD  capacitor  was  placed 
across  each  one  w'th  the  correct  damping  resistor  to  give  .7 
critical  damping.  A  UTC  low-pass  filter  L'-P-IO  was  placed  be¬ 
tween  each  of  the  seismometers  and  the  input  to  the  parametric 
amplifier.  This  action  insured  the  prevention  of  front-end 
ove  load  from  unwanted  high-frequency  data  prior  to  the  high- 
cut  filter  contained  in  the  later  stages  of  the  amplifier. 

A  seismometer  mass  centering  circuit  was  devised  utilizing 
the  calibrate  coil  in  each  of  the  sei  smom.jters  to  insure  the 


centering  of  the  mass  under  all  reasonable  positions  of  tha 
seismometer  and  with  motor  centering  voltages  as  low  as  6  V 
dc . 

It  was  found  that  of  the  six  seismometers  contained  in  the 
two  underwater  packages,  two  of  the  units  had  defective  sensing 
coils.  One  sensing  coil  was  shorted  to  the  calibration  coil 
and  the  other  sensing  coil  had  shorted  turns  reducing  the  coil 
resistance  from  40,000  ohms  to  28,000  ohms. 

Also  two  of  the  individual  centering  motors  had  to  be  re¬ 
placed  because  of  defective  gear  boxes.  This  condition  apparent 
ly  stemmed  from  the  external  gear  train  linkage  being  corpled 
too  tightly  to  the  mass  spring  system.  The  gear  teeth  were 
bottomed  into  each  other  causing  undue  friction  and  wear  on 
the  motor  gear  box  and  causing  the  motor  to  draw  more  current 
than  necessary  to  move  the  seismometer  centering  system.  This 
was  corrected  by  replacing  two  of  the  motors  and  placing  shims 
under  all  of  the  motor  mounts  as  well  as  slightly  moving  and 
lubricating  the  gear's. 

Modifications  of  the  Model  96  T^xas  Instrument's  tape  re¬ 
corders  were  also  made.  The  factory  recommended  changes  in¬ 
cluded  replacing  the  original  Barber-Col eman  BYQM-2953  6  V  dc 
tape  drive  motor  with  a  Sperry- Farragut  12  V  brushless  dc  motor 
Part  No.  490-229.  This  capstan  motor  change  necessitated  the 
design  and  installation  of  a  regulated  12  V  dc  power  supply  to 
maintain  accurate  tape  speed.  Also  included  among  the  recom¬ 
mended  changes,  was  rebuilding  the  motor  mounting  flywheel  plate 
and  setting  the  drive  linkage  wheel  bearing  in  an  epoxy  mounting 

Under  close  examinat'on  it  was  found  that  the  flywheels  of 
both  tape  recorders  were  cut  of  balance.  These  flywheels  must 
be  critically  balanced  becaused  under  dynamic  conditions  they 
rotate  at  1.6  RPS  If  left  uncorrected,  this  unbalance  could 
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wlthLW?'"’""  fron,  the  tape  recorder 

Within  the  desired  information  passband. 

for  th?,T"r"‘'  = 

the  Lake  Supertor  program  (10  to  13  days  per  drop)  it  was 

Lai  1  "  '",«'9ht.  The  amount  of  recording  tape  and  a 

n  incease  ,n  electronic  circuitry  also  Increased  the 
sphere  s  overall  weight.  The  Increased  weight  of  the  under- 
ater  package  together  with  the  problem  of  less  hr  yan cy 
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sphere  for  recovery  on  the  untethered  drop  This 
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system  response  Is  phase  matched  between  .5  and  8  cps. 

Augus7'97'7  ■'“"f- 

ugu  t  1967.  The  first  two  drops  are  to  be  tethered  drops  em- 

y  ng  one  ocean  bottom  seismograph  while  the  last  . -op  will 

u  1  ize  oth  seismographs,  one  being  a  free-fall  device  The 

latter  is  complicated  since  additional  batteries.  elect;on 

mponents.  and  reduced  buoyancy  (salt  to  fresh  water)  necL- 

a  es  adding  floatation  devices  and  a  different  anchor  ,e- 
lease  mechanism. 

Considerable  time  and  effort  was  devoted  to  this  task 
during  the  last  five  months.  The  first  unit  was  olaced  to 
operation  at  1943  EST  on  dune  6.  1967  with  the  asslstandof 
the  Coast  Guard  Cutter  Woodrush.  The  drop  point  was  in  the 
vicinity  Of  the  EARtV  RISE  shot  point.  The  unit  was  eco'  red 
0  oune  15.  1967.  The  total  recording  time  was  236  hou  and 
7  minntes.  Very  high  background  noise  was  encountered.  - 
'>-"ary  analyses  of  the  data  Indicates  an  average  peak  ground 
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displacement  of  apprcximate'’y  350  millimicrons  at  1  cps  and 
30  millimicrons  at  6  cps.  The  two  high  gain  data  channels 
for  each  component  were  overrecorded  most  of  the  time.  At 
three  times  during  the  recording  period  the  low  gain  channel 
was  overrecorded.  The  first  time  lasted  for  six  hours,  the 
'■•econd  for  60  hours  and  47  minutes  and  the  last  for  two  hours 
and  45  minutes.  Buildups  in  background  noise  of  5  minutes, 

1  minute,  and  25  minutes  preceded  these  periods.  Each  period 
was  followed  by  a  decay  in  background  noise  lasting  several 
hours  or  more.  At  this  time  it  is  undecided  whether  or  not 
this  is  due  to  natural  background  noise  or  instrumentation 
problems.  The  second  occurrence  began  so  abruptly  that  instru¬ 
mentation  problems  are  auspected.  Several  smaller  gradual 
buildup  and  decays  in  background  noise  appear  to  be  quite  legi¬ 
timate.  The  longest  one  started  at  0700  on  June  7,  1967 
(average  background  250  millimicrons),  reached  a  peak  at  approxi¬ 
mately  1630  on  June  8,  1967  (average  background  2520  milli¬ 
microns)  and  gradually  decayed  until  back  to  normal  at  2300 
on  June  9,  1967  (average  background  of  350  millimicrons).  Cur¬ 
rent  studies  are  being  conducted  to  correlate  this  information 
with  meteorological  data  and  with  the  land  seismic  station  that 
was  in  operation  during  the  same  time  period. 

At  least  eight  events  including  the  Buffalo  earthquake  of 
June  13,  1967  were  rec  rded  during  the  first  drop. 

A  second  drop  is  scheduled  for  July  5th  through  17th.  Modi¬ 
fications  in  the  amplifiers  and  damping  of  the  tape  recorder  base 
plate  mounting  system  will  be  conducted  before  this  drop. 
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2.6.  CRUSTAL  EFFECTS 


Preliminary  examination  of  the  EARLY  RISE  data  indicates 
a  considerable  amount  of  scattering  with  respect  to  normal 
attenuation.  This  phenomenon  has  also  been  characteristic  of 
explosion  data  previously  collected  in  the  same  geographical 
area  [2].  During  earlier  experiments  shot  size,  shot  location, 
and  station  location  were  all  varied  thus  making  it  difficult 
to  separate  shot  and  station  effects.  However,  with  the  EARLY 
RISE  project  both  shot  size  and  shot  location  were  fixed 
throughout  the  entire  experiment.  This  gives  one  an  excellent 
opportunity  to  examine  station  effects.  Our  attention  will  be 
given  primarily  to  the  analysis  of  body  waves  which  have  passed 
through  the  earth's  mantle. 

The  layers  of  the  earth's  crust  act  as  a  filter  with  re¬ 
spect  to  seismic  waves  arriving  at  a  given  station.  As  a  re¬ 
sult  the  motion  recorded  at  the  surface  depends  not  only  on 
the  frequency  x,ontent  of  the  explosive  source  but  also  on  the 
resDonse  of  the  layered  crust.  This  dependen  gives  us  a 
basis  for  a  method  of  determining  the  structure  of  the  crust 
under  each  station.  The  analysis  of  seismic  signals  for  this 
purpose  is  best  performed  in  the  frequency  domain.  In  order 
to  obtain  information  which  is  independent  of  the  frequency  of 
the  source,  the  spectrum  of  the  vertical  component  of  motion 
is  divided  by  the  spectrum  of  the  horizontal  component.  This 
ratio  represents  the  tangent  of  the  apparent  angle  of  emergence 
and  it  depends  only  on  the  angle  of  incidence  of  the  wave  and 
the  system  of  layers  below  the  recording  station.  This  is 
assuming,  however,  that  the  time  window  used  for  analysis  con¬ 
tains  only  reflectic.s  and  conversions  in  the  station  crust. 

Since  for  the  most  part  the  upper  mantle  is  homogeneous 
deep  focus  earthquakes  will  give  the  longest  time  window  for 
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the  direct  arrival.  In  this  case,  there  will  be  adequate  time 
for  the  most  significant  crustal  conversions  (to  Ci/)  and  reflec¬ 
tions  to  take  place  before  the  next  mantle  arrival.  With  body 
waves  from  an  explosive  source  recorded  at  less  than  1000  km 
the  situation  becomds  complicated  by  the  fact  that  there  will 
be  little  time  lag  between  the  mantle  and  the  first  crustal 
arrival . 

The  parameters  of  the  crustal  model  can  be  determined  by 
comparison  of  the  theoretical  and  observed  spectra  of  the  ratio 
of  the  components.  The  theoretical  curves  can  be  calculated 
easily  on  a  digital  oomputer  using  a  matrix  formulation  developed 
ty  Haskell  [3].  The  transfer  functions  for  both  the  horizontal 
and  the  vertical  component  can  be  calc  la  ted  as  a  function  of  fre¬ 
quency  given  a  model  where  the  compressional  velocity,  shear  velo¬ 
city,  density,  and  thickness  is  specified  for  each  layer  and  an 
angle  is  given  for  the  wave  incident  to  the  bottom  layer.  Their 
ratios  give  the  theoretical  tangent  of  the  apparent  angle  of 
emergence . 

Previous  investigations  [4],  [5],  and  [6]  have  been  limited 
to  using  lorig-period  body  waves  to  determine  a  one  or  two  layered 
crust.  The  dource  of  seismic  energy  has  been  limited  to  deep 
focus  earthquakes.  In  this  study  relatively  high  frequency  seis¬ 
mic  energy  (up  to  10.0  cps)  will  be  u^^e-^  to  determine  a  more 
complicated  crustal  model  with  special  emphasis  placed  upon 
near-surface  layering.  It  is  felt  that  the  upper  layering  con¬ 
tributes  most  to  the  nature  of  high  frequency  station  spectra. 

It  is  also  this  type  of  layering  that  could  vary  radically  along 
the  EARLY  RISE  station  line. 

In  order  to  use  seismic  data  for  crustal  determinations 
good  signal  to  noise  is  required  at  all  essential  frequencies. 

The  EARLY  RISE  data  has  adequate  signal  to  noise  as  far  as  the 
Mi  chi gan-Ohi 0  border.  After  this  point  the  acceptable  data  be¬ 
comes  spotty.  The  data  collected  in  the  Upper  Peninsula  of 

3fc 


I, 


! 


F- 

1 

1 

Michigan  had  exceptional  signal  to  noise  but  these  stations  lie 
too  close  to  the  shot  point  to  produce  the  necessary  separation 
of  prominent  phases.  Thus  the  greatest  portion  of  usable  data 
was  recorded  in  the  Lower  Peninsula.  The  geological  setting 
there  is  the  Michigan  Basin  structure  which  through  well  data 
is  known  in  great  detail.  One  Of  the  major  purposes  of  this 
study  is  to  try  to  correlate  the  major  station  anomalies  w  th 
the  Michigan  Basin  structure. 

Due  to  the  complexity  of  explosion  data  the  crustal  struc¬ 
ture  for  the  Botanical  Gardens  station  will  first  be  determined 
using  deep  focus  earthquakes.  These  should  ideally  be  100-600 
km  deep  and  at  a  distance  not  greater  than  6000  km.  A  search 
hds  been  made  of  the  U.S.  Coast  and  Geodetic  Survey  earthquake 
reports  for  this  type  of  event.  Several  have  been  fornd  and  a 
preliminary  frequency  analysis  is  being  made  of  each  Lo  ascer¬ 
tain  the  presence  and  distribution  of  high  frequency  energy. 

The  Brazi 1 i an-Peruvi an  earthquake  mentioned  earlier  appears  to 
be  acceptable  for  our  purpose.  We  hope  to  find  at  least  three 
or  four  such  earthquakes  with  some  high  frequency  content. 

These  will  then  te  digitized  and  analyzed  for  a  high 

resolution  of  frequency.  The  observed  tangent  of  the  apparent 
angle  of  emergence  will  then  be  compared  with  that  calculated 
from  a  proposed  crustal  model  for  the  Ann  Arbor  station. 

j  A  preliminary  model  is  suggested  on  the  basis  of  existing 

1  well  data.  The  model  consolidates  the  20  to  25  known  bosin 

layers  prior  to  the  basement  rock  ( PreCambr i an )  into  three  pri¬ 
mary  layers.  The  choice  was  mar  .  on  the  basis  of  the  degree  of 
velocity  contrast  between  adjoining  ’■'yers.  Depth  to  the  Pre 
Cambrian  structure  in  the  Ann  Arbor  area  is  approximately  5000 
I  feet  while  at  the  basin  center  it  is  over  13,000  feet. 
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PROPOSED  ANN  ARBOR  MODEL 


Layer 

Number 

Rock  Types 

a-Vel oci ty 
( km/sec ) 

8-Vel oci ty 
( km/sec ) 

Densi ty 
(qm/cmM 

Thi cknes  s 
(km) 

1 

Sandstone 
(porous ) 
Shale 

Some  Lime¬ 
stone 

2.8 

1.6 

2.4 

0.29 

2 

Dol omi tes  , 
Salts 

5.5 

3.2 

2.8 

0.70 

3 

Limestone , 
Sandstone 
(hard ) 

3.9 

2,3 

2.5 

0.76 

4 

PreCambri an 
(Grani te ) 

6.2 

3.6 

2.6 

10.0 

5 

Unknown 

7.7 

4.^- 

2.7 

20.0 

6 

(mantl e) 

8.3 

4.  « 

3.0 

-  - 

Calculations  using  the  Thomson-Haskel 1  matrix  formulation 
and  the  above  model  indicate  that  the  basin  layers  would  be 
sensitive  to  high  frequency  seismic  energy.  Assuming  that  a 
consistent  model  can  be  found  for  the  well  station  if]  attempt 
will  be  made  to  determine  the  crustal  profile  along  the  EARLY 
RISE  station  line  in  the  Lower  Peninsula.  EARLY  RISE  shot 
spectrum  will  all  br  normalized  using  the  Botanical  Gardens 
records.  Any  complexities  of  the  explosion  data  will  then  be 
determined.  The  variations  in  the  ratio  spectra  from  that 
predicted  by  the  earthquake  determined  model  will  be  due  chief¬ 
ly  to  body  waves  passing  through  the  lower  crust.  If  any  dis¬ 
crepancies  can  be  accounted  for,  corrections  can  be  made  to  the 
shot  spectrum  for  each  station  and  the  individual  station  crusts 
can  then  be  fitted. 
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Some  difficulty  may  occur  in  trying  to  determine  a  basin 
model.  Fernandez  and  others  have  pointed  out  the  complexities 
that  arise  when  transfer  functions  are  computed  at  high  ;re- 
quency.  This  is  also  compounded  with  the  difficulty  in  trying 
to  determine  a  rather  complex  model.  It  is  expected  at  least 
that  this  endeavor  will  give  some  knowledge  to  the  relation¬ 
ship  between  the  amount  of  data,  its  frequency  content,  and 
the  degree  of  crustal  complexity  which  can  reliably  be  re¬ 
solved  using  the  spectrum  of  body  waves. 
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Appendix  I 

CHRONOLOGICAL  ANNOTATFD  BIBLIOGRAPHY  OF  REPORTS 
PREPARED  UNDER  CONTRACT  AF  49(63C)-1759  an 
AND  IN  PART  BY  AF  49(638)-n70 


David  E.  Willis  and  John  M.  DeNoyer,  "Seismic  Attenuation  and 
Spectral  Measurements  from  the  Lake  Superior  Experiments,"  Am . 
Geophys.  Union  Monograph  #10,  "The  Earth  Beneath  the  Continents," 
1966. 

Detailed  attenuation  studies  of  the  seismic  data  obtained 
from  the  Lake  Superior  experiment  showed  frequent  rapid  changes 
in  the  spectra  and  maximum  amplitudes.  In  attempting  to  analyze 
these  data,  a  technique  was  developed  that  permits  the  presen¬ 
tation  of  the  seismic  data  in  three  dimensions:  distance,  fre¬ 
quency  content,  and  particle  velocity  amplitude.  Contoured  iso¬ 
particle-velocity  maps  are  thus  prepared  which  show  that  many 
of  the  rapid  changes  in  spectra  and  amplitude  are  a  part  of 
significant  trends  extending  from  more  than  one  shot  point. 

The  results  show  that  conditions  in  the  Source  region  for  under¬ 
water  shots  strongly  influence  the  character  of  the  seismic  waves 
recorded  at  considerable  distandes.  Residual  iso-particle- 
velocity  maps  can  help  to  isolate  the  shot  anomalies  and  provide 
considerable  insight  to  the  overall  pattern  of  the  attenuation  of 
the  seismic  waves.  Using  an  equation  of  the  type  A  =  AQR-nE-ctfR 
to  express  the  attenuation  of  the  seismic  waves,  values'^of 
a  -  3  to  5  X  10‘'’sec/km  and  n  =  1  or  3/2  were  formed  to  give  rea¬ 
sonable  fits  to  the  observed  first  compressional  wave  arrivals. 
Longer  range  data  obtained  from  the  October  1964  Lake  Superior 
shots  indicated  a  value  of  a  =  2  x  10"^  sec/km  and  n  =  3/2. 

John  M.  DeNoyer,  G.  E.  Frantti,  and  D.  E.  Willis,  "Short  Note  on 
Underwater  Sound  Measurements  from  the  Lake  Superior  Experiment," 
Am.  Geophys.  Union  Monograph  #10,  "The  Earth  Beneath  the  Cont’- 
nents,"  1^66. - 

Underwater  sound  recordings  of  the  Lake  Superior  shots  showed 
that  the  predominant  energy  was  contained  in  a  band  between  60 
and  250  cps.  Late  arrivals  were  found  that  correlate  with  re¬ 
flections  froin  the  main  shoreline,  numerous  islands,  and  a  few 
prominent  sho-'ls.  By  recording  on  magnetic  tape  and  by  using 
various  passband  filters  on  the  laboratory  playback,  it  was 
found  that  more  precise  times  could  be  determined  for  the  water 
wave  arrivals. 


Willis,  "A  Note  on  The  Comparison  of  Short-Period  Seis- 
for  Recording  Teleseismic  Events,"  Earthquake  Notes, 

No.  2,  June  1966. 
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Six  short-period  seismometers  of  different  makes  anc  modeli, 
which  are  in  common  use,  were  used  to  record  several  underground 
nuclear  detonations  at  teleseismic  distances.  A  comparison  of 
the  seismogrems  obtained  showed  no  clear  superiority  of  one  type 
seismometer  over  anotner.  However,  arrival  times  for  a  given 
phase  recorded  by  the  various  s.:i smrmeters  showed  discrepancies 
which  could  introduce  errors  in  crustal  and  upper  mantle  struc¬ 
tural  determinations. 

P.  L.  Jackson,  "Correlation  Function  Spatial  Filtering  with  In¬ 
coherent  Light,"  Applied  Optics,  July  1967. 

The  method  of  lensless  correl ography  (Meyer-Eppler  eff->ct) 
has  been  extended  to  image  one  of  the  transparencies  which  are 
being  correlated.  An  optical  filter  is  placed  in  the  correla¬ 
tion  plane,  so  that  the  image  of  the  transparency  represents 
the  "time-domain"  analog  of  a  complex  frequency  filter. 

P.  L.  Jackson,  "Truncations  and  Phase  Relationships  of  Sinusoids," 
J.  Geophys.  Res.,  Vol .  72,  No.  4,  1967,  pp,  1400-1403. 

The  theoretical  basis  of  frequency  shifts  in  the  Fourier 
transform  of  short  segments  is  given,  ihese  frequency  shifts 
arise  from  constructive  and  destructive  interferences  between 
the  positive  and  negative  frequency  components  the  Fourier 
tran5fo<"m.  Surprisingly,  similar  shifts  are  found  in  electronic 
bandpass  filtering. 

D.  E.  Willis,  "A  Long-Range  Reversed  Sei.smic  Profile  in  the 
Eastern  United  St^tp' ,"  abstract  for  presentation  at  the  I.U.G.G. 
Meeting,  Zurich,  196/. 

A  series  of  underwater  high-explosive  shots  fired  during  the 
past  five  years  have  provided  data  for  a  reversed  seismic  pro¬ 
file  approximately  2000  km  in  length.  These  shots  include  the 
1962  offshore  North  Carolina  experi'^'^nt ,  the  1965  ECOOE  series, 
the  CHASE  III,  IV,  afid  VI!  events,  ai.J  the  Lake  Superior  shots 
of  1953,  1964,  and  1966.  Close  station  intervals  were  obtained 
along  the  profile.  Travel-time  curves,  frequency  analyses  and 
attenuation  measurements  are  presented  for  these  recordings.  An 
intermediate  crustal  layer  with  a  velocity  of  7.7  km/sec  was  dis¬ 
closed  on  the  northwestern  portion  of  the  p/ofile.  Higher  P 
velocities  and  a  thicker  crust  was  also  indicated.  Smaller  " 
travel-time  residuals  across  the  Michigan  Basin  were  observed  on 
both  profiles.  The  travel-time  data  indicated  a  possible  moun¬ 
tain  root  system  under  the  Appalachians.  Contoaued  particle- 
velocity  maps  made  from  the  spectral  data  displosed  amplitude 
anomalies  that  could  be  correlated  with  conditions  in  the  source 
region.  These  results  aid  in  the  interpolation  of  t!ie  attenua¬ 
tion  data. 
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